The QCD string tension curve, the ferromagnetic magnetization, 
and the quark-antiquark confining potential at finite Temperature 
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We study the string tension as a function of temperature, fitting the SU(3) lattice QCD finite 
temperature free energy potentials computed by the Bielefeld group. We compare the string tension 
points with order parameter curves of ferromagnets, superconductors or string models, all related to 
confinement. We also compare the SU(3) string tension with the one of SU(2) Lattice QCD. With the 
curve providing the best fit to the finite temperature string tensions, the spontaneous magnetization 
curve, we then show how to include finite temperature, in the state of the art confining and chiral 
invariant quark models. 
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I. INTRODUCTION 

Motivated by the discovery of the Quark-Gluon Plasma 
at the CMS at CERN, and RHIC at BNL laboratories 
and by the future experiments at RHIC-II, LHC and 
FAIR, quark model computations are presently address- 
ing with more and more detail chiral symmetry breaking 
and deconfincmcnt at finite temperature or finite den- 
sity. Here wc show how to upgrade the confining quark- 
antiquark potential, with a string tension a function of 
the temperature T. 

In the 80's a new quark model, the Chiral Quark 
Model (xQM), was developed including not only 

the quark confining potential of the early quark models, 
but also the spontaneous breaking of chiral symmetry 
of the Nambu and Jona-Lasinio model. Although the 
xQM can only be approximately derived from QCD in 
the Coulomb gauge hamiltonian formalism, although it 
has not yet been fully calibrated to the finest details of 
the hadronic spectrum, and although Lorentz invariance 
is only present in the kinetic energy, it is nevertheless the 
only QCD inspired model able to explicitly include both 
quark confinement and quark-antiquark vacuum conden- 
sation. Thus the xQM is the most comprehensive QCD 
model, adequate to study any possible system of quarks, 
gluons, or hadrons, and to provide, at least, a qualitative 
answer to the questions one may ask to QCD. 

So far the xQ-M studies of chiral symmetry breaking 
at finite temperature or density have assumed a tempera- 
ture independent string tension, i. e. the same confining 
potential at all temperatures, since the work of Yaouanc, 
Oliver, Pene , Raynal, Jarfi and Lazrak Q. Within this 
approximation, many studies of chiral symmetry have 
found that chiral symmetry breaking is maintained at 
all temperatures. The first calculations at finite temper- 
ature or at finite density are already 20 years old 0, Q , 
but presently finite temperature or finite density are in- 
teresting many authors, including the Sao Paulo group 
[1, 0] , the Graz group [1] , the Pittsburgh group Q and 
others. Only Lo and Swanson using ring diagrams, i. 
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FIG. f: We show examples finite temperature static quark- 
antiquark potentials, in particular the T < Tc and T > Tc 
Lattice QCD data for the free energy Fi, thanks to [lol - 
[l3]01af Kaczmarek et al. The solid line represents the T — 
static quark-antiquark potential. In this paper we discuss the 
use of the free energy as a finite temperature quark-antiquark 
potential. 



e. quark loops beyond the ladder approximation, has 
been able to find chiral symmetry restoration at a finite 
critical temperature Tc, since the quark loops effectively 
affect the confining potential. 
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FIG. 2: Fits of the Bielefeld's group free energy aX T < Tc, 
here in the case of T = 0.90Tc, in the top we use a linear fit 
and in the bottom we use a linear + coulomb fit. We cut the 
low distance, or infrared, part of the lattice data shown in 
Fig. [T] in such a way that the fit is independent as possible 
on the cutoff. Similar cuts and fits are applied for the other 
temperatures. 



Nevertheless the finite temperature static quark- 
antiquark free and internal energies at some finite tem- 
peratures T have already been computed in SU(3) lat- 
tice QCD, by the Bielefeld group |10l - [l4l |. and they arc 
available to be included in the xQM. In particular, the 
first part of the quark-antiquark potential to sufi'er finite 
temperature effects is the long distance confining part. 
In what concerns the short range Coulomb part of the 
quark-antiquark potential, it survives up to any temper- 
ature. Moreover the short-range Coulomb potential is 
included in the ultraviolet rcnormalization program of 
QCD, thus it is relevant for the quark mass rcnormal- 
ization, but can be neglected for chiral symmetry break- 
ing [l^. While the Coulomb potential at finite tempera- 
ture is quite relevant for the heavy quarkonia, since even 
slightly above the deconfinement critical temperature Tc 
the J/ ij} and rjc remain bound [iB, II3l , for the binding 
of light quarks the confining part of the quark-antiquark 
potential, parametrized by the linear string tension, is 
crucial. 

Moreover, although the Polyakov loop is presently the 
preferred order parameter for deconfinement [Tsl . [l9j . 
other order parameters may be used, (20l . [2l| and in par- 



T/Ta 




IR cutoff xV dof 


0.00 


1.545 ± 0.002 


1.0 




0.90 


0.861± 0.072 


3.6 


0.52 


0.94 


0.587 ± 0.028 


4.3 


1.96 


0.98 


0.429 ± 0.005 


4.1 


0.22 


1.00 


0.00 ± 0.000 







TABLE I: Fits of the string tension in units of Tc" V(0) ^''^ 
for the longer distance part of the Kackzmarek et al lattice 
data, using the ansatz F{r) = c + err. 



ticular the string tension is a possible order parameter 
for the deconfinement phase transition. 

Thus here we specialize in the finite temperature string 
22| . and show how to include the results of fi- 



tension 

nite T Lattice QCD in the xQ^i including finite T in 
the string tension of the confining part of the quark- 
antiquark potential. This works continues the study of 
the SU(3) string tension, computed for the first time by 
the Bielefeld Group, [13 • 

In Section II we fit the SU(3) string tension from the 
long distance part of the recent static quark-antiquark 
free energy from the lattice QCD results of the Bielefeld 
Group, and we compare the SU(3) string tension with 
the SU(2) string tension. In section III we compare the 
string tensions with the order parameter curves of ferro- 
magnets, superconductors or string models, all related to 
confinement, and find that the curve closer to the SU(3) 
string tension is the magnetization curve. In Section IV 
we use the curve inspired in the magnetization of fer- 
romagnets, to also compute the entropy string tension, 
the internal energy string tension, and to estimate the 
finite T string tension of the quark-antiquark potential. 
In Section V, we conclude. 



II. FITTING THE SU(3) AND SU(2) LATTICE 
QCD FINITE T STRING TENSIONS WITH 
DIFFERENT ANSATZE 

This works continues the study of the SU(3) string ten- 
sion, computed for the first time by the Bielefeld Group, 
[2^ . in a seminal paper also discussing in detail the fi- 
nite temperature phenomenology. In this 2000 paper, the 
Bielefeld group computes the string tension below Tc uti- 
lizing colour averaged (over the colour singlet and colour 
octet) Polyakov loops, to extract colour averaged free en- 
ergies. The results then obtained confirmed a 1st order 
phase transition for the deconfinement phase transition, 
as expected for SU(3) since the string tension was found 
to remain finite at T < Tc, while it is vanishing at T > Tc, 
with a small discontinuity at T = Tc- 

More recently, the Bielefeld group computed free ener- 
gies for the colour singlet quark-antiquark static system, 
utilizing gauge fixing. Here we fit the string tension of 
these free energies, [l(]| - |l4| recently computed by Doring, 
Hubner, Kaczmarek, Karsch, Vogt and Zantow in SU(3) 
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TABLE II: Fits of the string tension in units of T^" V(0)~^'^^ 
for the longer distance part of the Kackzmarek et al lattice 
data, using the ansatz F{r) = c — ctj^ + ar. 



Lattice QCD by , for T < T^, sec Fig. [H For T > Tc 
the fit is trivial since the confinement potential, and the 
string tension, vanish. To extract the string tension as a 
function of the temperature cr(T) we try different ansatze 
to fit the free energy, and we cut the low distance part 
with an infrared cutoff. We choose the lowest possible 
cutoff, in the region where our fit is stable for changes of 
the cutoff. We also check that our fit is not too far from 
the of x/dof ~ 1. As different ansatze to fit the long 
distance part of the potential, we use a constant and a 
linear term, 

F{r) =c{T) + <j{T)r (1) 
or a constant, a coulomb and a linear terms 

F{r) ^ c{T) - ^ + a{T)r (2) 
r 

or, as in the SU(3) Bielefeld fit [12] and in the similar 
SU(2) fit of reference [2^, a constant, a linear and a 
Logarithmic terms 



F{r) = c(T) - a{T) \og{r) + a{T)r 



(3) 



and in all of them we check the stability of the string 
tension. We apply these ansatze to the free energies at 
the different temperatures, and this is illustrated if Fig. 
[5]for the case of T = 0.90Tc. It occurs that the fit with a 
Logarithmic term of Eq. ([3]) is not sufficiently stable, in 
the sense that whenever we increase the infrared cutoff, 
the string tension cr change, and thus we abandon this 
ansatz. The other two fits of Eqs. ([T]) and ([2]) are stable, 
in the sense that a distance rcut can be used as an in- 
frared cutoff, and that if we increase this infrared cutoff 
distance, the string tension is essentially unchanged. 

The details on the fits are show in Tables |T] and |lll 
The simpler constant plus linear fit of of Eq. ([T]) needs 
a higher cutoff rcut than the constant plus linear plus 
Coulomb fit of eq. ([3]). The two different sets of string 
tensions are depicted in Fig. [3l Notice that to the finite 
temperatures of T = 0.907^, T = OMT^, and T = 0.98Tc 
we add the cases of T = with the string tension cr(0) 
and oi T > Tc with a vanishing string tension cr(T') = 0. 

It occurs however that the string tensions computed 
with the two different ansatze of Eqs. ([T]) and ^ dif- 
fer, even when they are normalized to Notice that 



0.8 



0.4 



0.2 



0.6 0.8 



1.2 Tc 



1.0* 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 



1.0 



1.2 Tc 



FIG. 3: The results of the fits, present in Table [T] (top) and 
in Table [n] (bottom) , can be summarized in a dimensionless 
curve of a{T)/a{0) as a function of T/Tc- The bottom fit, of 
a constant plus a Coulomb plus a linear term has a greater 
stability that the simpler top fit, and smaller error bars. 



the r = potential has the Luscher term and the 
same term is expected to also exist at finite temperature, 
the ansatz including a Coulomb term of Eq. ([2]) has a 
broader stability for changes of the infrared cutoff rcut, 
and produces points with a smoother alignment. Thus 
of all our fits of the string tension of the Bielefeld free 
energies, the ansatz of Eq. ©, with the fits of Table HIl 
and Fig. [3] (bottom) appears to be the best. 

In Fig. Uwe plot our fits, present in TableUand in Ta- 
ble |lll together with the string tensions a that were fitted 
by the Bielefeld Group in 2000 They are all consis- 
tent, with our linear fit slightly larger than the Bielefeld 
fits in the year 2000, and with our linear plus Coulomb fit 
only slightly lower than the Bielefeld fit. In Fig. [5] we also 
compare the SU(3) string tension critical curve (bullets) 
with the SU( 2) string tension critical curve (squares) of 
reference [2J]. The SU(2) sting tension is smaller than 
the SU(3) one, and it corresponds to a larger critical ex- 
ponent than the critical exponent of SU(3). 
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FIG. 4: The results of the fits, present in Table U and in 
Table nil together with the SU(3) fits in the year 2000 of the 
Bielefeld Group [l^] , summarized in a dimensionless curve of 
cr(T)/(j(0) as a function of T/Tc. 



FIG. 5: Comparing the SU(3) string tension critical curve 
(bullets, squares and diamonds) of Fig. |4l with the SU(2) 
string tension critical curve (triangles) of Ref . [23| . The criti- 
cal curve of SU(2) is flatter, and it corresponds to a larger crit- 
ical exponent 0.63 than the critical exponent ~ 0.5 of SU(3). 



III. FERROMAGNETIC, SUPERCONDUCTOR 
AND STRING INSPIRED ANSATZE FOR THE 
QCD STRING TENSION CURVE 

The string tension a is crucial for the spontaneous 
breaking of chiral symmetry, and for the light hadron 
spectra. Thus it is important to know the string tension 
cr(T) at all temperatures, staring from T = and not 
just at temperatures of the order of Tc. Since in Section 
II the string tension was only computed for a few tem- 
peratures, close to Tc, we now compare the string tension 
to similar other curves or order parameters, with the aim 
to propose an ansatz for cr(T). 

We remark that we expect for QCD a 1st order phase 
transition with a discontinuity at T = Tc, but with only 
a weak discontinuity. Thus we may use as ansatze for the 
string tension curve, in the case where we are interested 
in all the interval T S [0,Tc] and not just in a narrow 
neighbourhood of Tc, second order parameter curves. In- 
deed the points in Fig. [3] are close to second order, or 
continuous phase transition, with a critical exponent in 
the order parameter. Notice that we do not claim here 
that the transition is not the expected 1st order one, we 
are simply concerned with the fit of the potential to be 
used in the quark model. 

To find an ansatz for the string tension curve we first 
study order parameter curves of physical systems related 
to confinement. The Ising model not only is a model 
of ferromagnetism, but it is also a model of confinement, 
and in particular the critical exponent of the SU(2) string 
tension is similar to the three-dimensional Ising model 
exponent. Another model of confinement is the dual su- 
perconductor model, since the colour electric flux tubes 



existing between confined quarks and antiquarks are con- 
fined as vortices in superconductors. And colour electric 
flux tube conflnement can be approximated, in the limit 
of thin flux tubes, by string models. Thus we may inspire 
ourselves in ferromagnetic materials, in the Ising model, 
in superconductors either in the BCS model or in the 
Ginzburg-Landau model, or in string models, to suggest 
ansatze for the string tension curve. 

We order the different ansatze used in these related 
systems, from the simpler to define, to the more sophis- 
ticated one. The coherence length of a Ginsburg-Landau 
superconductor [^J] is, 

the string tension of a finite T string [1^, extrapolated 
from T ~ Tc to T ~ 0, is, 




the empirical fit to the penetration length of a supercon- 
ductor [iBi is. 




the spontaneous magnetization of a fcrromagnct [27l | is 
the solution of the algebraic equation, 

^^tanhf^JL), (7) 
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FIG. 6: We show the string tension curves —. — " ,,,, as a func- 
tion of Y- of the parametric ansatze of Eq. (|4]) in a dashed 
hne, Eq. ([SJ in a dot-dashed hne , and Eq. ((SJ in a dotted 
Una. The solution of Eq. ((8]) for a large Debye frequency wd 
is similar to the dashed line. 
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FIG. 7: The critical curve for as a function of ^ is show 

1^1 sat J- c 

in a dotted solid line (one dot per each of our numerical fixed 
point solution). We find the critical curve solving Eq. ((T)) 
with the fixed point method. The solution is essentially a flat 
constant curve for low T but for T ~ Tc it behaves like the 
square root as in eq. also depicted here. 



and the mass gap of a BCS superconductor [1^ is a so- 
lution of the integral equation, 



.giV(O) 



dw 



tanh ■ 



2f 



(8) 



We compare the curves of Eqs. (jH), ([5|), © in Fig. [51 
and we show how to plot the curve of Eq. ([7]) in Fig. [T] 
The mass gap solution of Eq. (|5]) is not unique since it 
depends on two parameters, the density of states gN{0) 
and the Dcbyc frequency wd- In the limit of a very small 
Wd, Eq. ([8]) is equivalent to Eq. ([7]). The opposite limit 
Wd » gN{0) occurs in real type I superconductors j26j . 
in which case the curve A(T), in dimensionless units, is 
quite similar to the circular curve of Eq. ([5]). Thus we 
do not plot the solution of Eq. ([8|) in a separate curve. 

In Fig. m we plot the Magnetization, solution of Eq. 
(O, together with the string tensions a that fitted from 
SU(3) Lattice QCD data. The solution of Eq. © pro- 
vides our best fit of the string SU(3) tensions. We also 
show in Fig. [S] our second best ansatz, the one of the 
empirical penetration length of a ferromagnet in eq. ^ . 
In Fig. iniwe also compare the magnetization curve with 
the SU(3) string tension critical curve computed in Rcf. 
m, and with the fit 



ct(O) 



1.21i 



0.990 



(9) 



used in Ref. [22| to measure the finite string tension at 
as an evidence for a 1st order phase transition. Both the 
fit of Eq. © , similar to Eq. ^ but with a different norm 



and with a small temperature shift of lO^^Tc, showing 
evidence for the 1st order phase transition, and the SU(3) 
string are quite close to our magnetization curve for T 
close to Tc- But the magnetization curve is the only one 
that extends up to the correct solution at T = 0. The 
magnetization curve is as close to our string tension curve 
as it is close to the magnetization of a real ferromagnet. 

Since the magnetization curve fits so well the SU(3) 
string tensions, we now briefiy review how it is derived. 
The ferromagnetic spontaneous magnetization curve as a 
function of temperature is a text book curve, detailed for 
instance in the Feynman Lectures on Physics 27[. Let 



us briefly review the simplest case, with spins in the ^ 
representation. Then, in a magnetic field B, each spin 
has only two possible states with energies, ±^o B where 
/iQ = f ^ is not the chemical potential but it is the 
magnetic moment of a quantum spin. This leads to the 
magnetization. 



M = NfiQ tanh 



kT 



(10) 



The case that interests is the case of a ferromagnet, when 
there is no external field, and the magnetic mean field 
affecting the spin is proportional to the magnetization 
M . NfiQ is the saturation magnetization Msat- We then 
get the algebraic Eq. ([7]) equation for the magnetization 
curve. 



M 



tanh 



Tc M 

T Msal 



the Eq. providing the best ansatz for our string tension 
curve. 
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FIG. 8: Comparing the ferromagnet magnetization M/Msat 
critical curve (solid line) of Fig. [7] with the SU(3) string ten- 
sion CT/cr(0) (bullets, squares and diamonds) of Fig. [21 both 
as a function of T ~ Tc. In this plot, M/Msat is quite close 
to ct/(j(0). We also show (dotted line) the empirical curve of 
Eq. (|6} which provides the second best fit of all the ones we 
tried. 
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FIG. 9: Comparing the magnetization curve (solid line) 
with the SU(3) string tension critical curve (bullets) com- 
puted in Ref. [S^]. We also show (dashed line) the fit 
1.21(1 - 0.990(T/Tc)^)^/^ used in Ref. [H] to measure the 
finite string tension at Tc as an evidence for a 1st order phase 
transition. 



To analyze the curve solution to Eq. ([7]), we start by 
getting the extreme points of the curve solution to eq. 
([7]) , noticing that tanhoo = 1, thus when T = we get 
= 1. Also when the only solution oi x = tanhx is 

X = thus when T = Tc we have jj^ = as expected. 
Moreover, expanding the hyperbolic tangent close to T ~ 
Tc, we get , 



M ^ Tc M 1 fTc M V 




This shows that in this case the critical exponent for 
the magnetization M is i. This corresponds to a sec- 
ond phase transition since the magnetization is the first 
derivative of the free energy with regards to the chem- 
ical potential. For a complete solution of Eq. ([7]), we 
use the fixed point method, with 10^ iterations for each 
temperature T . In Fig. [7] we show the solution of Eq. 
([7]) obtained with the fixed point expansion, and also the 
approximate solution shown in of Eq. (jlip for T ~ Tc- 

It is striking that the curve -^^^^^p-^, which is derived 
from a simple two-state quantum spin system, not only 
fits the magnetization curve of many fcrromagncts, but 
also fits essentially the string tension curve of finite tem- 
perature quenched SU(3) QCD. Possibly this occurs since 
confinement can also be partially understood in the sim- 
plest two-state Ising Model. In particular Digal, Fortu- 
nato and Pctreczky [2^ computed the string tension of 
quenched SU(2), also with fits of the free energy, and 



found that the SU(2) string tension has the same critical 
exponent 0.63 of the 3-dimensional Ising model. Thus in 
Fig. [5] we also compare the SU(3) string tension criti- 
cal curve (bullets) with the SU(2) string tension critical 
curve (squares) of reference |23| . The critical curve of 
SU(2) is lower than the SU(3) curve, and it corresponds 
to a slightly larger critical exponent than the critical ex- 
ponent ~ 0.5 of SU(3), obtained when neglecting the first 
order discontinuity at T = Tc. Nevertheless the similar- 
ity between our SU(3) string tension and the the SU(2) 
and Ising model string tensions suggests that the SU(3) 
confinement is relatively simple. Indeed Lattice QCD is 
able to simulate comprehensively hadronic physics using 
a relatively small number of 100 to 1000 configurations, 
and this is only possible if SU(3) confinement has rela- 
tively few degrees of freedom. It may also be relevant for 
the understanding of SU(3) confinement, that the mean 
field model of ferromagnetism provides a better fit to the 
SU(3) string tension curve than the Ising spin-spin inter- 
action model. Also, when comparing with a the circular 
curve of Eq. ([5]), the SU(3) string tension curve is fiatter 
at the origin, and thus closer to the curve of Eq. ([6]), 
while the SU(2) string tension curve is steeper and thus 
closer to the curve of Eq. 

Thus, based on our numerical results, and also confi- 
dent in the relative simplicity of SU(3) QCD confinement, 
and on the small discontinuity of the string tension a{T) 
at Tc we may use, in what concerns finite temperature 
quark models, as an ansatz for the string tension, 

a{T/Tc)^a{Q)^WEA . (12) 
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FIG. 10: A free illustration of a static flux tube, where we 
compare the force / to the pressure P, and the length r to 
the volume V . 



FIG. 11: Free energy (solid line) compared with the inter- 
nal energy (dashed line) (left), entropy divergent at T = Tc 
(right). 



IV. THE FREE ENERGY, THE INTERNAL 
ENERGY, AND THE FINITE TEMPERATURE 
QUARK-ANTIQUARK CONFINING 
POTENTIAL 

In the xQM, the potential energy V is used in the mass 
gap equation to compute the quark constituent mass and 
chiral symmetry breaking, and in the Salpeter equation 
to compute the hadron spectrum. However, in Lattice 
QCD, first the quark-antiquark free energy F is com- 
puted from the pair of Polyakov loops P, 



Tr{P(0)P^(r)} ^ Ccxp 



-F 



(13) 



and it is the free energy of the a static quark-antiquark 
pair that we have studied in Sections II and III. The 
potential energy of the quark-antiquark pair is related to 
the quark-antiquark force. 



dV 



-fdr 



(14) 



where a linear flux tube is assumed as in Fig. 1101 and 
we can study the thermodynamics of the flux tube using 
the equivalence f dr = p dV for the work of the quark- 
antiquark pair. Then we get the usual thermodynamic 
relations, relevant at finite T, for the free energy F, 



F{r) = -fdrSdT , 
and for the internal energy E, 

E{r) ^-fdr + TdS 



(15) 



(16) 



The quark potential V is identical to the free energy F 
in a reversible isothermal transformation with dT = , 



and is identical to the internal energy E' in a reversible 
adiabatic transformation with TdS = 0. 

Since we have a good ansatz for the confining part of 
F, i. e. for the string tension, determined in Section III, 
we may compute the string tension of the entropy S, 



S 



dF 



(17) 



and then also compute the string tension of the internal 
energy, 



E^ F + TS 



(18) 




In Fig. [TTl we compare the different string tensions. No- 
tice that the entropy diverges when T — >■ Tc, and this is 
well known from the computations of the Bielefeld Group 
at temperatures close to Tc- In particular the divergence 
is, 



(19) 



and thus the internal energy E also diverges when T — > 
Tc, this is clear in Fig. [TT] 

Now, in a hadron, or in the vacuum, it is not to- 
tally clear what energy, either F or i? to use as a the 
quark-antiquark potential. A main assumption of the 
xQM, and of any quark model, is that the static quark- 
antiquark computed in Lattice QCD can be also be used, 
at least qualitatively, for light quarks. One assumes that 
the quark moves relatively slowly compared to the flux 
tube formation typical scale, and that the gluonic flux 
tube adapts to the quark and antiquark positions. This 
is an adiabatic assumption, in the quantum mechanical 
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FIG. 12: Flattening the potential combining some internal 
energy to the free energy (left), flattening the free (solid line) 
and internal (dashed line) energies with a square root (right). 



sense and not in the thermodynamic sense. Thus with 
the mechanical adiabatic assumption of a slow motion 
of the flux tube, the quark movement is sufficiently slow 
in the heat bath of the hot medium at temperature T, 
and we can assume that the flux tube transformation is 
isothermal. Then the potential energy is well approxi- 
mated by the free energy V ^ F. Nevertheless, possi- 
bly the quarks move too fast for a completely isothermal 
transformation of the flux tube, and then the flux tube 
transformation may have a small contribution from the 
internal energy E, identical to the potential energy when 
no heat is exchanged with the heat bath, in that case. 



V = {l-uj)F + LoE , 



(20) 



where w should be a small number. In that case the 
potential dependence in the temperature T is flatter for 
T < Tc- In particular we also depict in Fig. [12] the cases 
where 0.1 < to < 0.3 where the string tension is nearly 
flat up to close to T — Tc- 

Notice however that there is an evidence contrary to 
the use of the internal energy E for light quarks. The 
problem with the internal energy is not only it's diver- 
gence at T = Tc, but also that it is larger than the free 
energy, E > F. Yamamoto, Suganuma and lida found 
in Lattice QCD that a the presence of light quark does 
reduces the string tension [23|. Thus, presently, the best 
model for the string tension at finite T remains the one 
of Eqs. (0 and (HH). 



V. CONCLUSION 

We fit the finite temperature string tension from the 
free energies of a quark-antiquark pair cornputed in SU(3) 
lattice QCD by the Bielefeld group |10l4l4j . We find a 
good ansatz for the string tension at all temperatures, the 
magnetization curve of a ferromagnet. While evidence for 
a 1st order phase transition exists [22|, the difference to 
the 2nd order magnetization curve is small, and therefore 
the magnetization curve is a good ansatz for the quark- 
antiquark potential at finite temperature, adequate to be 
included in the xQM. 

We find that a constant string tension may also be 
arguably used at T < Tc- Thus at T < Tc the quark 
model computations present in the literature are accept- 
able. The main arguments are that the free energy F 
string tension is already relatively flat up to T close to 
Tc and that the light hadron spectrum scales with ^/a 
which is flatter than the string tension, as in Fig. [T^ 
Moreover, if a contribution from the internal energy E 
string tension can be justified, then the string tension is 
further flattened up to T close to Tc- And the Irst or- 
der discontinuity, although small, also flattens the string 
tension. 

However, at T > Tc the string tension vanishes, and in 
this sense the xQM calculations present in the literature 
must be corrected. To encompass all temperatures, the 
free energy string tension of Eqs. ([7]) and ([T^ should 
be used in the xQM, or at least a step function with a 
transition at Tc could be approximately used. 

Nevertheless, the relative instability of the flts of the 
finite temperature string tension, and the importance of 
the finite T Coulomb or logarithmic potentials suggest 
that eventually the finite T xQM will have to go beyond 
the simple use of a long range linear potential, including 
also medium range and short range potentials [H, [1^ . 

In what concerns the Lattice QCD studies of the free 
energy at finite temperature, it would be interesting if 
more results would be further computed. In particular 
the subtle 1st order SU(3) discontinuity at T = Tc [2^ . 
deserves more detailed studies in the region of T ~ Tc. 
Moreover at smaller temperatures, T s]0,0.8Tc[ there is 
no Lattice QCD data. Thus we anticipate that the re- 
search of finite T quark-antiquark potentials will remain 
very interesting in the future. 
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